WEYL-TITCHMARSH M-FUNCTION ASYMPTOTICS FOR 
MATRIX- VALUED SCHRODINGER OPERATORS 



STEVE CLARK AND FRITZ GESZTESY 



Abstract. We explicitly determine the high-energy asymptotics for Weyl- 
Titchmarsh matrices associated with general matrix-valued Schrodinger oper- 
ators on a half-line. 



1. Introduction 

The high-energy asymptotics z — > oo of scalar-valued Weyl-Titchmarsh func- 
tions m+(z,Xo) associated with general half-line Schrodinger differential expres- 
sions of the type —-^2 + q{x), with real-valued coefficients q^ n ' £ i 1 ([x ,c]) for 
some n £ No(= N U {0}) and all c > xq, received enormous attention over the 
past three decades as can be inferred, for instance from [6], [7], [8], [9], [11], [25], 
[31], [32], [33], [46]-[49], [50], [51], [61], [62], [86] and the literature therein. Hence 
it may perhaps come as a surprise that the corresponding matrix extension of 
this problem, considering general matrix-valued differential expressions of the type 
—Im^ri +Q(x), with I m the identity matrix in C m , m £ N, and Q(x) a self-adjoint 
matrix satisfying Q(") £ L 1 ^, c]) roxm for some n £ No and all c > Xq, received 
hardly any attention at all. In fact, we are only aware of a total of two papers 
devoted to Weyl-Titchmarsh M-function asymptotics for general half-line matrix 
Schrodinger operators: one by Clark [20] in the case m = 2, n — 0, and a second 
one, an apparently unpublished manuscript by Atkinson [7] in the case m £ N, 
n = 0. Both authors are solely concerned with determining the leading term in 
the asymptotic high-energy expansion of the Weyl-Titchmarsh matrix M+(z,xo) 
as z — > oo without investigating higher-order expansion coefficients. (It should be 
noted that this observation discounts papers in the special scattering theoretic case 
concerned with short-range coefficients QS n > £ L 1 ([x ,oo); (1 + |a;|)c?x) mxm , where 
straightforward iterations of Volterra-type integral equations yield the asymptotic 
high-energy expansion of M + (z, xq) to any order, cf. Lemma 4.1.) A quick look at 
Atkinson's preprint reveals the nontrivial nature of determining even the leading- 
order term of the asymptotic expansion of M + (z,xo) as z — > oo, 

M+(z,xo) = iz 1 l 2 I m + o(\z\ 1 ' 2 ) (1.1) 

2 — >OQ 

for z in sectors in the open upper half-plane due to a variety of variable transforma- 
tions, which finally reduce the problem to the study of a managable Riccati system 
of differential equations. 

Our principal motivation in studying this problem stems from our general inter- 
est in operator- valued Herglotz functions (cf. [12], [37], [38], [39], [40], [41], [44], [85]) 
and its possible applications in the areas of inverse spectral theory and completely 
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integrable systems. More precisely, using higher-order expansions of the type (1.1), 
one can prove trace formulas for Q(x) and certain higher-order differential polyno- 
mials in Q(x), following an approach pioneered in [42] (see also [35], [36]). These 
trace formulas, in turn, then can be used to prove various results in inverse spectral 
theory for matrix- valued Schrodinger operators H — —I m 4-^ + Q in L 2 (R)"\ For 
instance, using the principal result of this paper, Theorem 4.6, and its straightfor- 
ward application to the asymptotic high-energy expansion of the diagonal Green's 
matrix G(z, x, x) — (H — z)^ 1 (x, x) of H , the following matrix-valued extension of 
a classical uniqueness result of Borg [10] has been obtained in [22]. 

Theorem 1.1 ([22]). Suppose H is reflectionless (e.g., Q is periodic and H has 
uniform spectral multiplicity 2m) and spectrum [E , oo) for some E a £ M. Then 

Q(x) = E I m for a.e. x £ M. (1.2) 

For related results see, for instance, [2], [13], [14], [15], [17], [18], [19], [23], 
[26], [59], [67], [68], [69], [77], [82]. Incidentally, the higher-order differential poly- 
nomials in Q(x) just alluded to represent the Korteweg-deVries (KdV) invariants 
(i.e., densities associated with KdV conservation laws) and hence open the link to 
infinite-dimensional completely integrable systems (cf. [3], [16], [27], [28], [29], [30], 
[34], [45], [70], [71], [72], [75], [83], [84] and the references therein). 

In Section 2 we briefly recall some of the basic notions of Weyl-Titchmarsh theory 
for singular matrix Schrodinger operators on a half-line as developed in detail by 
Hinton and Shaw [54]-[58] (see also [5], [50], [60], [63], [64], [65], [66], [76], [78], [79], 
[80], [83], [89]). In fact, most of these references deal with more general singular 
Hamiltonian systems and hence we here specialize this material to the matrix- 
valued Schrodinger operator case at hand. The subsequent section is devoted to 
a proof of the leading-order asymptotic high-energy expansion (1.1) of M + (z, xq), 
originally due to Atkinson. Since the result is a fundamental ingredient for our 
principal Section 4, and to the best of our knowledge, his manuscript [7] remained 
unpublished, we provide a detailed description of his approach in Section 3. Finally, 
Section 4 develops a systematic higher-order high-energy asymptotic expansion of 
M+(z,xo), combining Atkinson's result in Section 3 with matrix-valued extensions 
of some methods based on an associated Riccati equation. More precisely, following 
a technique in [43] in the scalar- valued context, we show how to derive the general 
high-energy asymptotic expansion of M + (z,xo) as z — » oo by combining Atkin- 
son's leading-order term in (1.1) and the corresponding asymptotic expansion of 
M_|_(z, xq) in the special case where Q has compact support. 

Analogous results for Dirac-type operators are in preparation [21]. 

2. Weyl-Titchmarsh Matrices 

In this section we briefly recall the Weyl-Titchmarsh theory for matrix-valued 
Schrodinger operators on half- lines (x , oo) for some xq £ M. 

Throughout this note all m x m matrices M £ C mxm will be considered over 
the field of complex numbers C. Moreover, I p denotes the identity matrix in C pxp 
for p £ N, M* the adjoint (i.e., complex conjugate transpose), M* the transpose of 
the matrix M, and AC([o, b]) denotes the set of absolutely continuous functions on 
[a,b\. 

The basic assumption of this paper will be the following. 
Hypothesis 2.1. Let m £ N and suppose Q = Q* £ L 1 ([xo,c]) mxm for all c > x . 
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Given Hypothesis 2.1 we consider the matrix- valued Schrodinger operator H + 
mL 2 ({x Q ,oo)) m , 

d 2 

dom(H + ) = {g G L 2 ((x ,oo)) m \g,g' G AC{[x 0l c]) m for all c > a; ; lim #(a;) 

s.-a.b.c. at oo; (-/ m .g" + Qg) G L 2 ((x , oo)) m }, 

where the abbreviation s.-a.b.c. denotes a fixed self-adjoint boundary condition 
at oo throughout this paper if the differential expression — I, 



(2.1) 
0. 



Q{x) is not 

in the limit point case at oo (cf. the paragraph following Definition 2.7) and is 
disregarded in the case where —I m ^i + Q{x) is in the limit point case at oo. 
Naturally associated with the operator H + is the matrix Schrodinger equation 

— ip"(z, x) + Q(x)tp(z, x) = zip(z, x), z^C,x>x , (2.2a) 

where z represents a spectral parameter and tp(z,x) is assumed to satisfy 

ip'(z, •) e AC([x ,c]) mxr for all c> 0, (2.2b) 

with the value of r depending upon the context involved. Equation (2.2a) may also 
be expressed as the system 

.M (z, x)' = (zA(x) + B{x))^{z 1 x), z G C, x > x , 



where 



and where 



W(z,x) 
A(x) = 



ip(z,x) 
tp'{z,x) 







-In 









B(x) 



-Q{x) 




(2.3a) 
(2.3b) 
(2.3c) 

(2.3d) 



G AC([x Q ,c}) 2mxr for all c> 0; 

again, with the value of r depending upon the context involved. 

Next we turn to Weyl-Titchmarsh theory associated with (2.1) and briefly recall 
some of the results developed by Hinton and Shaw in a series of papers devoted 
to spectral theory of (singular) Hamiltonian systems [54]-[58] (see also [5], [50], 
[60], [63], [64], [65], [66], [76], [78], [79], [80], [83], [89]). While they discuss much 
more general systems, we here confine ourselves to the special case of matrix- valued 
Schrodinger operators governed by Hypothesis 2.1. 

Let *&(z, x, xo) be a normalized fundamental matrix of solutions of (2.3) at some 
xo G R; that is, ^(z, x, xa) is of the type 

$>(z,x,x ) = [i>j,k{z,x,x )]\ k=1 

= [@{z,x,x ) $(z,x,x )] 
*(z,x ,a;o) = I 2m , 

where Q(z, x, xq), $>(z, x, xq) G C 2mxm , and where 9(z,x,xo) and (j>(z, x,xq) rep- 
resent a fundamental system ofmxm matrix-valued solutions of (2.2), entire with 
respect to z G C, and normalized according to (2.4b); that is, 



6(z, x, x ) 
9'(z,x,x ) 



(f>(z,x,x ) 
f>'(z,x,x ) 



(2.4a) 
(2.4b) 



9(z,x ,x ) = cf)'(z,xo,x ) 



6'(z,x Q ,x ) = <p(z,x ,x ) = 0. (2.4c) 
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Let = [Pi 2 }; with U 2 e C mxm . We assume 

ker (/?*) fl ker (/3 2 ) = {0}, or equivalently that rank(/3) = m, (2.5a) 
and that either 

±{l/2i)0J0* = ±Im(/J 2 /J*) > 0, or Im^/?*) = 0. (2.5b) 

Here we denote, as usual, Im(M) = (M - M*)/(2i) and Re(M) = (M + M*)/2. 
We now prove the following result concerning $(z, x, Xo). 

Lemma 2.2. Let <&(z,x,x ) be defined in (2.4), and let G C mx2m saiis/y (2.5). 
If 0$>(z, c, xo) = Pi<f>(z, c, xo) + 02<fi'( z , c, xo) is singular for c > x and Im(/?2/9*) ^ 
0, i/ierc Im(z) «s 0. // 0&(z, c, x ) is singular for c > x and Im(/3 2 /3i) = 0, £/ien 
Im(z) = 0. 

Proof. Let <f(z,x) satisfy (2.3). Then, 

&(z,c)*M(z,c) -&(z,x )*J&(z,x ) = 2ilm(z) / ds<f (z, s)M(s)f (z, s). (2.6a) 

Let \P(z,x) = &(z,x,Xo) and note that 

V(z,xo)*M(z,x ) = [0 I m ]J[0 ImY = 0. (2.6b) 

Moreover, note that 

&(z, c)*M{z, c) = <f>'(z, c, x )*<f>(z, c, x ) - <f>(z, c, x )*<j>'(z, c, x ). (2.6c) 

Now suppose that 0<S>(z,c,x o )v = for v G C m , u ^ 0. If Im(/? 2 /?*) ^ 0, then 
(j>'(z,c,Xo)v = —02 1 0i4>{z,c,xq)v. Thus by (2.6) we obtain, 

v*(f>(z, c, xo)*/3 2 " 1 [Im(/3 2 /3J')]/32~ 1 0(2, c, x )v 

= — Im(z) / ds v*cp(z, s, x o )*0(z, s, x )w, (2.7) 

J Xo 

by which we obtain the first part of the lemma. However, if = Im(/3 2 /3 1 ) = 
(l/2i)0J0* then by (2.5a) there is a w G C m such that 

$(z,c,x )u = J(3*w. (2.8) 

For this case, by (2.6) we obtain 

= Im(z) / ds v*<t>(z, s, x )*(f>(z, s, x )v, (2.9) 

J x 

which implies that Im(z) = 0. □ 
One can also prove the following result. 

Lemma 2.3. Let Q(z, x, x ) and x, x ) be defined in (2.4), and let G C mx2m 
satisfy (2.5). Fore > x , /3$(z,c, x ) is singular if and only if z is an eigenvalue 
for the regular boundary value problem given by (2.3) together with the boundary 
conditions 

[I m 0}&{z, x ) = 0, 0$(z, c) = 0. (2.10) 

Lemmas 2.2 and 2.3 show that with appropriate values of z e C one may define 
a certain m x m matrix. 
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Definition 2.4. Let Q(z,x) and $(z,x) be denned in (2.4), and let (3 G C mx2m 
satisfy (2.5). Then, for c > x and for (3<&{z,c,Xo) nonsingular, let 

M(z,c,x ,(3) = -(f3^(z,c 1 x Q ))- 1 (l3e(z,c,x )). (2.11) 

M(z, c, xo, (3) is said to be the Weyl-Titchmarsh M-function for the regular bound- 
ary value problem described in Lemma 2.3. Also, let D(z,c,Xo) denote the set of 
all M(z, c, xo, (3) defined in (2.11) together with the assumption (2.5b). 

About M(z, c, xo, (3), we prove the following result. 

Lemma 2.5. Let c > Xo- If Im(/?2/3*) > 0, then M(z,c,xo, (3) is well-defined if 
Im(z) > 0; i/Im(/?2/3*) < 0, then M(z, c, xo, (3) is well- defined for Im(2:) < 0; 
and z/Im(/32/9i) = then M(z, c, Xo, [3) is well-defined for z G C\R. Moreover, if 
Im(/3 2 /3*) | and Im(z) ^ 7 then 

Im(±M(z,c,x ,f3))>0. (2.12) 

Proof. Except for the conditions given in (2.12), the claims made in the statement 
of this lemma are an immediate consequence of Lemma 2.2. 
To verify the conditions given in (2.12), we begin by letting 



U(z, x, x ) = [®{z, x, x ) ${z, x, x )] 



M 



(2.13) 



where Q(z, x, x ) and $(z, x, x ) are defined in (2.4), and where M G C mxm . Now, 
let ^(z,x) — U(z,x,xo) in (2.6a) and note that 

*(z, x )*.m(z, x a ) = 2Tm(M). (2.14) 

Moreover, if M = M(z, c, xo, (3), then also note that 

f3V(z,c) = (3U{z,c,xo) = 0. (2.15) 

If (l/2i)(3J/3* = Im(/3 2 /3*) = 0, then by (2.5a) and (2.15) there is a matrix 
C G C mxm such that V(z,c) = U{z,c,x ) = J[3*C. Thus, 

*(z,c)*J*(z,c) =C*[3J[3*C = 0, (2.16) 

and (2.12) follows immediately from (2.14), (2.16) and (2.6a) when Im(z) ^ 0. On 
the other hand, if Im(/3 2 /3i) ^ 0, then (2.15) implies that u'(z,c) = — P^ 1 [3iu(z , c) . 
As a result, 

*(z,c)*J*(z,c) = -2iu{z,c)*^ 1 luY{l32l3l)l3*^ 1 u{z 1 c), (2.17) 
and again (2.12) follows from (2.14), (2.17), and (2.6a) when Im(z) ^0. □ 

Remark 2.6. By Lemma 2.5, M(z, c, xo, 13) is a Herglotz function of rank m. More- 
over, D(z, c, xo) is contained in the standard Weyl disk of matrices associated with 
(2.2) or (2.3) for fixed values of z, c, and xo (cf. [21]). These disks nest with respect 
to increasing values of c; that is, D(z, Ci, Xo) Q D(z, c\, Xo) whenever xo < c\ < ci. 
This leads to the existence of a closed limiting set D + (z, xo), which we call a lim- 
iting disk and which corresponds to the standard limiting Weyl disk (cf. [21], [54], 
[56], [57], [76]). 

Definition 2.7. Let D + (z, x ) denote the limit, as c — > oo, of the nested collection 
of sets D(z,c,xo) given in Definition 2.4, and let M + {z,xo) denote a point of the 
limit disk D + (z,xq) associated with (2.2). 
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We note that D + (z,x ) consists of a unique matrix if and only if —I m d 2 /dx 2 + 
Q(x) is in the limit point case at +00 and that D + (z, xo) consists of a closed, convex, 
and bounded set with nonempty interior if and only if —I m d 2 /dx 2 + Q(x) is in the 
limit circle case at +00. In both of these cases, we note that elements of dD + (z, x ) 
represent half-line Weyl-Titchmarsh matrices where each such element is associated 
with the construction of self-adjoint operator H + on the half-line. However, for 
those operators that fall between the limit point and limit circle cases, Hinton and 
Schneider have noted that not every element of dD + (z, xq) is a half- line Weyl- 
Titchmarsh matrix and have characterized those elements of the boundary that are 
(cf. [52], [53]). 

For convenience of the reader we quickly summarize the principal results on 
half-line Weyl-Titchmarsh matrices in the following theorem. 

Theorem 2.8 ([4], [12], [44], [54], [55], [58], [64]). Assume Hypothesis 2.1, x G 
R, z G C\M, and let M + (z,x ) be a half-line Weyl-Titchmarsh matrix. Then 

(i) M + (z,x ) is a matrix-valued Herglotz function of maximal rank. In particular, 

Im(M + (z,a;o)) > 0, z G C+, (2.18) 

M + (z,x ) = M + (z,x y, (2.19) 

rank(M+(z, x )) = m, (2.20) 

limM + (A + ie,x ) exists for a.e. A G M. (2.21) 

Isolated poles of M + (z,x ) and — M + (z, x^)^ 1 are at most of first order, are real, 
and have a nonpositive residue. 

(ii) M + (z,xq) admits the representations 

M + {z,x )=F + (x )+ f dQ + (X,x )((X-z)- 1 -A(l + A 2 )-!) (2.22) 

JR 

= cxp(c+{x ) + J dXE + (X,x )((\ - z)- 1 - A(l + A 2 )- 1 )^), (2.23) 



where 



F + (x )=F + (x y, / J'^jyU oo, (2.24) 
C+(x ) = C+(x y, <- + (■, x ) <I m &.e. (2.25) 



Moreover, 



fi+((A, n],x ) = limlim- / dvluv{M + (v + ie,x Q )), (2.26) 

<5i0 slO 7T J x+S 

E+(X,x ) = limilm(ln(M + (A + ie,xo)) for a.e. X G R. (2.27) 

ej.0 7T 

(in) With <j)(z,x,xo) and 9(z,x,x ) given in (2.4), define the m x m matrices 

u + (z,x,x ) = 0(z,x,x o ) + (f>(z,x,xo)M + (z,x ), x, x n G M, z G C\M. (2.28) 
Then 

u +{ z i ■) x o) £ L 2 {(xq, oo)) mxm , (2.29) 
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u + (z,x,x n ) is invertible, u + (z,x,x ) satisfies the boundary condition of H + at 
infinity (if any), and 

/>oo 

Im(M + (z,x )) —lm(z) / dx u + {z, x, x )*u + (z, x, x ). (2.30) 

3. Atkinson's Leading-Order Argument 

The purpose of this section is to prove the following result due to Atkinson. 

Theorem 3.1. (Atkinson [7].) Assume Hypothesis 2.1 and denote by C e C C+ the 
open sector with vertex at zero, symmetry axis along the positive imaginary axis, 
and opening angle e, with < e < n/2. Let M + (z,xo) be either the unique limit 
point or a point of the limit disk D + (z, x ) associated with (2.2). Then 

M+(z,x ) = iz 1/2 I m + o(\z\ 1/2 ). (3.1) 

\z\ — >oo 
z£C e 

More precisely, Atkinson proved the analog of (3.1) in [7] for more general Sturm- 
Liouville-type situations rather than matrix-valued Schrodingcr operators. To the 
best of our knowledge, Atkinson's manuscript [7] remains unpublished. Since the 
result is crucial for our principal Section 4, we will specialize his arguments to the 
simplified situation at hand and provide a complete account of his approach (at 
times going beyond some of the arguments sketched in [7] ) . 

One of our tools in studying M + (z, x n ) will be to relate it to a matrix Riccati-type 
equation with solution M + (z, x) defined as follows. We first note that 

M + (z,x ) = M+(z,.T )u+(z,a;o) _1 , z G C\R, (3.2) 

where u + (z, ■ ) G L 2 ([x , oo)) mxm is a nonnormalized solution of (2.2) satisfying 

u+(z, x) = u + (z, x, x n )C + , (3.3) 

with u + defined in (2.28), and with C + G C mxm a nonsingular matrix. Varying 
the reference point xq we may define 

M + (z,x) = u' + (z,a;)M + (z,a;) _:L , z G C\R, x > x Q . (3.4) 

Then M + (z,x) is independent of the chosen normalization C+ and is well-known 
(see, e.g., [22], [35], [60]) to satisfy (3.5) below. 

Lemma 3.2. Assume Hypothesis 2.1, suppose z G C\M, x G [x ,oo), and define 
M + (z,x) as in (3.4). Then M + (z,x) satisfies the standard Riccati-type equation, 

M' + (z,x)+M+(z,x) 2 =Q(x)- zl m . (3.5) 

Proof. Differentiating (3.4) with respect to x and inserting u + (z,x) + (Q(x) — 
z)u + {z,x) = immediately yields (3.5). □ 

Utilizing 

*(z, x, x\) = *(z, x, x a )^(z, x Q , xi) and *(z, x , x 1 )^f(z, x x ,x ) = I 2m , (3.6) 
(cf. (2.4)) one computes for all x\ G M 

9(z, x, x\) + <p(z, x, xi)M + (z, x\) 

= u + (z,x,x o )[0(z,x o ,xi) + <j)(z,x ,xi)M + (z,xi)}, z G C\R. (3.7) 
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Since u+(z,;x ) G L 2 ([R, oo)) mxm for all R G M, the left-hand side of (3.7) is 
in L 2 ([R, oo)) mxm . Thus, for M + (z,x\) defined according to (3.4), one concludes 
that 



M + (z,x\) G D + (z,xi) for all iri G 



(3i 



since M + (z, x ) € Z?+(z,a;o) by hypothesis. This justifies our choice of notation 
M + (z, x) when compared with the notation employed in Definition 2.7. 

Atkinson's proof of Theorem 3.1, which we follow very closely, is first con- 
cerned with the relationship between D(z,c,xq), described in Definition 2.4, and 
D n (z, c, xo), a Riccati disk, which we now define. 

Definition 3.3. With Q(x) given in Hypothesis 2.1 and z G C+, D n (z,c,x ) is 
defined to be the set of all M(z,x ) G C mxm where Im(M(z, x )) > 0, and such 
that the solution of 



&(z,x) = ^[l m + 0(z,x) I m -0(Z,X) 



—l\Z 



. 1 / 2 {zI m -Q{ X )) 

i\z\ l / 2 I, 

I m + <d{z,x) 
I m - $(z,x) 



(3.9a) 



ti(z,x ) = (I m +i\z\- 1 / 2 M{z,x ))(I m -i\z\- 1 l 2 M{z,xv))- 1 

= (z|z| 1/2 / m - M(z,x ))(i\z\^ 2 I m + M(z,x Q ))-\ (3.9b) 



satisfies 



d{z, x)*d(z 1 x) < I m for all x G [xq,c]. 



(3.10) 



Note that there is a correspondence between elements of D n (z, c, Xq) and certain 
solutions of (3.9a). If -d(z,x) satisfies (3.9a) and (3.10) then note that a matrix 
M(z, xo) G D n (z, c, xo) is defined by the initial data d(z, xq) if $(z, Xo)*$(z, x ) < 
I m . Note also that (3.9a) can be written as 





i\z\ 1/2 i„ 



{zA + B) 



-i\A- 1/2 Ir, 









(3.11) 



Remark 3.4. We emphasize that (3.9a) formally results from the Cayley-type trans- 
formation, 



t?(z, x) = {i\z\ x / 2 I m - M(z, x))(i\z\^ 2 I m + M{z, x))-\ 



(3.12) 



where M(z,x) satisfies the Riccati-type equation (3.5). In the scalar context this 
corresponds to a conformal mapping of the complex upper half-plane to the unit 
disk. Moreover, we note that Atkinson proves the asymptotic result of Theorem 3.1 
for elements of D n (z,c,xo). Recall that D + (z,xq) C D(z,c, xq) (cf. Remark 2.6). 
Theorem 3.1 then follows upon showing that D(z, c, xo) Q D K (z, c, x ) for c > x - 
This containment is shown in Theorem 3.6. 

Lemma 3.5. Ifu(z,x) is defined by (2.13) with M G D(z,c,Xo), then u(z,x) — 
i\z\^ 1 / 2 u'{z,x) is invertible for all z G C +) and all x G [x ,c]. 
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= A/3* - foft = [ft ft] 



(3.14) 



Proof. With M G D(z,c,x ) and with U(z,c) defined in (2.13), U(z,c) satisfies 
(2.15). If Im(/J 2 /?i) > 0, then u'(z,c) = -f5^ X {i\u{z, c). U(z,c) has rank to, hence 
it (2, c) is nonsingular. As a result, 

lm(u(z, c)*u'(z, c)) = «(«, c^/^Im^/^K" 1 ^ c) > 0. (3.13) 
However, if Im(/3 2 /3*) = 0, then 

-fill 

By (2.5a), rank[/3i /3 2 ] = rank[/3 2 - /3i]* = m; thus, by (2.15) and (3.14) there is a 
w e C mxm such that £/(z, c) = [f3 2 - /3i]*w. From this, we obtain 

Im(u(z, c)*u'(z, c)) - w*Im(/? 2/ 3*)?« = 0. (3.15) 

Hence for /3 which satisfy (2.5), 

lm{u(z,c)*u'(z,c))>Q. (3.16) 

For all x G [x , c] and z G C+, note that 

(Im(u(z, x) V(z, a;)))' = -lm(z)u(z, x)*u{z, x) < 0. (3.17) 

Together, (3.16) and (3.17) imply that 

Im(u(z, x)*u'(z, a;)) > for all a; G [a;o, c], z G C+. (3.18) 

If for some £ G [xo,c] there is an / G C m , / 7^ 0, such that u(z,£)f = 
i\z\- 1 / 2 u'(z^)f, then 

-i^ 1 / 2 /*"^ 0/ = /*«(*, 0*«'(*, 0/- (3-19) 

Together, (3.18) and (3.19) imply that f*u(z, g)*u(z, £)/ < 0, and hence that 
0/ = 0/ = 0- Given the uniqueness of solutions for system (2.3), we con- 
clude that = u(z, x )f = /; thereby producing a contradiction which completes 
the proof. □ 

Theorem 3.6. D(z, c, xo) C D n (z, c, xo) /or a// c > xo and z G C+. 

Proof. With z G C + and M — M(z,c,x o ,0) G D(z,c,x Q ), we see that Im(M) > 
by (2.12). With u(z,x) defined in (2.13), then by Lemma 3.5, we can define 

d{z,x) = (u(z,x)+i\z\- 1/2 u'(z,x))(u{z,x)-i\z\- 1/2 u'(z,x))- 1 (3.20) 

for x G [xo, c]. As defined, fl(z,x) satisfies (3.9b). Then for x G [x c], (3.18) 
implies that 

2i|z| _1 / 2 {w(z, x)*u'(z, x) - u*'(z, x)u(z, x)} 
= ~4\z\- 1/2 lm(u*(z,x)u'(z,x)) < 0. ( ' > '~" 
This is equivalent to 

(u(z, x)* - i\z\- 1/2 u'{ Zl x)*)(u(z, x) + i\z\- 1/2 u'(z, x)) 
< {u{z,x)* +i\z\- 1 / 2 u'{z,x)*){u{z,x) -i\z\-^ 2 u'(z,x)). 

Given the invertibility of u(z,x) — i|z|~ 1 / 2 u'(z, x) shown in Lemma 3.5, we infer 
that z9(z,x) satisfies (3.10). The proof of this theorem will be completed upon 
showing that $(z,x) satisfies (3.9a). 
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First, observe that 

tf(z,x)(u{z,x) - i\z\~ 1/2 u'(z,x)) = u(z,x) + i\z\- 1/2 u'(z,x), (3.23) 
and hence that 

I m + 1? = 2u(u - ilzl" 1 / 2 ^)" 1 , I m -# = -2i\z\-^ 2 u'(u - i^- 1 ' 2 ^)- 1 . (3.24) 
Upon differentiating (3.23) and using the fact that u(z,x) satisfies (2.2), we obtain 
0'(u - i|*r V V) - (I m - !?)«' - i\z\- x l 2 {I m + d){zl m - Q)u. (3.25) 
By (3.24), it follows that 

= "2^p72 (7m ~ ~ ^2~~ (7m + ??)(z/m ~ 0)(/m + ^' (3 - 26) 
which is precisely (3.9a). □ 

The following relation holds for the sets in Theorem 3.6, that is, 

D(z, c, x ) = D n (z, c, xq) for all c > x and z e C+. (3.27) 

This will be discussed in detail in [21]. 

An associated system for (3.9) is obtained by introducing a change of independent 
variable: With z 6 C e , c 6 [x , oo), and x E [xo, c], let 

i = (x-x )|z| 1/2 . (3.28) 

Introducing 

<ys(z,i) = i?(z,x), (3.29) 



with x and t related as in (3.28), (3.9) becomes 

1 

2 



</>'(-M) = 7J [Im + <P(z,t) I m -<P(z,t)] 



i\z\-\zI m -Q{t)) 
Urn 



I m + ip(z, t) 
I m - <p(z, t) 



(3.30a) 



<p(z, 0) = (i\z\V 2 I m - M{z, x ))(i\z\^ 2 I m + M(z, xo))" 1 , (3.30b) 
and (3.10) becomes 

<p(z,t)*<p(z,t) <I m for all* <= [0, (c - x )|z| 1/2 ]. (3.31) 
Note that in (3.30a) 

Q(t)=Q(x + t\z\- 1 / 2 ). (3.32) 

In (3.30) and (3.31), we have a set of conditions equivalent to (3.9) and (3.10) for 
D n (z 7 c, xo) given in Dcfintion 3.3. 

Now consider a sequence, z„ £ C El such that \z n \ — > oo as n — > oo and such that 

< s < 6 n = arg (z n ) <ir~e. (3.33) 

By choosing an appropriate subsequence, we may assume that 

tf„-»5e [e,7r-e]. (3.34) 

Let <p(z n ,t) denote a corresponding sequence of functions that satisfy (3.30a) and 
(3.31), with initial data, tp(z n ,0), defined by (3.30b) for a sequence of points 
M(z n ,Xo), neN, where each M(z n ,Xo) is chosen to lie in the disk D n (z, l7 c, x ). 
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Note that as z n — ► oo, the intervals described in (3.31) eventually cover all compact 
subintervals of M.+ . Given the uniform boundedness of tp n {t) = <p(z n ,t) described 
in (3.31), we assume, upon passing to an appropriate subsequence still denoted by 
f n (0), that 

9?n(0) = <p{z n ,0) — ► ipo(S), as n — ► oo, (3.35) 
and as a consequence, that 

¥>o(<*)Vo(<J) < Jm- (3.36) 

With ipo(5) defined in (3.35) as \z n \ — > oo, we consider a limiting system associ- 
ated with (3.30): 



V '(8,t)=-[l m +v(S,t) I m -r)(5,t)] 

t > 0, 



-ie iS I„ 






il m 



T](5,t) 

n(S,t) 



(3.37a) 
(3.37b) 

(3.38) 



Theorem 3.7. The solution f](6,t) of (3.37) satisfies 

ri(S,t)*ri(5,t) < I m 

for < t < oo. Moreover, the solutions <p n (t) = <p(z n ,t) of (3.30) converge to 
r](S,t) uniformly on [0, T] for every T > 0, as n — > oo. 

Proof. Let T e M+ be the greatest value such that (3.38) holds for t e [0, T]. Wc 
show that (3.38) must hold for some [0, T'] with T' > T, thus proving T = oo. 

The solution of (3.37), rj(5, t), presumed to be defined on [0, T], can be continued 
onto some [0, T'\ with T' > T; r](5,t) then satisfies 

r](5,t)*r](6,t) < k 2 I m (3.39) 

for < t < T and for some k > 1. 

For brevity, let ip' n = G n ((p n ,t) denote (3.30a) with z — z n , and let rj' = H(n,t) 
denote (3.37a) in the following. Integrating (3.37) and (3.30), we obtain 

<Pn(t) - V(fi,t) = ip n (0) - ipo(S) + / ds{G n (n,s) - H(rj,s)} 

Jo 

+ / ds{G n {ip n ,s) - G n (r), s)}. 



(3.40) 



Note that 



G n ( V , s) - H( V , s) = -i(e lS - e iS ")(I m + V {5, s)) 2 



+ -i\z n \- x (I m + n(8, s))Q(s){I m + n(6, *)), 



and by (3.39), as n — > oo and for t E [0, T'], that 



\i{e iS - e iS ") f ds (I m +i 1 {5, s )) 2 = 0{5-5 n ), 



(3.41) 



(3.42) 



and by (3.28), (3.32), and (3.39), that 



ds {I m +r](S,s))Q{s){I m +r]{6,s)) 
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= \z n \- 1/2 0ij dx\\Q(x)\\) (3.43) 

= o(|z„p 1/2 ) as n -> oo. (3.44) 
Together with (3.35), we observe that as n — > oo (3.40) is equivalent to 

<p n {t) - T](S, t) = o(l) + / ds(G n {<p n , s) - G n {r), a)) (3.45) 
Jo 

uniformly so, for < t < T . Since \\ip n \\ < 1 and |?i|| < k for t e [0, T'}, 

\\G n {<p n ,8)-G n { V , 8 )\\ < ^(3 + k)\\ l p n ( s )- v (S, S )\\(2+\z n \- 1 \\Q(s)\\). (3.46) 

Thus (3.40) yields 

\\<p n (t)-ri(5,t)\\ 

<o(l)+j o ds 1(3 + A)|| V „(*)- V(S,8)\\(2 + \z n \- 1 \\Q{s)\\). 

In light of (3.44), an application of Gronwall's inequality to (3.47) yields 

(p n (t) - Tj(5, £) — ► as 7i — > oo (3.48) 

uniformly for < t < T". For n sufficiently large and t G [0, T'], (p n (t) satisfies 
(3.31) with z = z n , and r](S,t) satisfies (3.39), hence r](S,t) satisfies (3.38) for 
t G [0, T'\, where V >T. □ 

The following computation identifies fo{S) and proves that r](S,t) is constant 
with respect to t > 0. 



(3.47) 



Corollary 3.8. 



V (5, t) = V (5, o) = MS) = \ - cM fJfl ^ *>o. (3.49) 

1 + cxp(io/2) 

Proof. Utilizing the standard connection between the explicit exponential solutions 
of the second-order Schrodinger equation (2.2a) and the Riccati equation (3.5) (in 
the special case Q(x) — 0, cf. (3.2) and Lemma 3.2), performing a conformal map 
of the type (3.12), and the variable transformations (3.28) and (3.29) then yields 
the following solution for (3.30a), 

<p(z,t) =( - (i\z\ 1/2 - iz 1/2 ){M (z) +iz 1/2 I m ) +cxp(-2itexp(i(5/2)) x 

x (i^/ 2 + lz 1 / 2 )(M (z) - iz 1 ' 2 ^)) x 

x ( - (i\z\ 1/2 + iz 1/2 ){M (z) + iz 1/2 I m ) + exp(-2itexp(i6/2)) x 

x (ilzl 1 / 2 - iz^ 2 )(M (z) - iz^Im))' 1 , (3.50) 

associated with the general initial condition 

<p(z, 0) = (i\z\ l ' 2 I m - Mo(z)) (i\z\ l ' 2 I m + M (z)y l (3.51) 

for some M (z) G C mxm with Im(M (z)) > 0, z G C+. Since by hypothesis 
< 5 < 7r, the exponential terms in (3.50) enforce 

\\<p(z,t)\\ > 1 as t t oo (3.52) 

unless 

M Q {z) = iz 1 ' 2 I m , (3.53) 
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implying (3.49). □ 

Given these facts we proceed to the 
Conclusion of the proof of Theorem 3.1. 

With M(z n ,Xo) E D n (z,c,x ), and for S as in (3.34) let 

_ 1 - cxp(^/2) 
C[S) l + exp(iJ/2)- (3 ' 54) 
From (3.30b), (3. 35), (3.48), and (3.49), wc infer that 

(I m +i\z n \- 1 / 2 M(z ni x Q ))(I m -i\z n \- 1 ' 2 M{z n ,x ))- 1 -» C(S)I m as n -> oo. 

(3.55) 

As a consequence, 

•i 1-1/2 w/ \ C(5) - 1 + o(l) r C(<f)-1,., „,, T 

,|*.| ^ M (z w , g0 ) = + l + Q \^ m= e^fTi (1 + 0(l))/ ro 

= -^/ 2 (l + o(l))/ m . (3.56) 
Thus, by (3.54) and (3.56) we conclude that 

M(z n ,x )=iz 1 J 2 (l + o(l))I mi (3.57) 

and hence by Theorem 3.6, that (3.1) holds. □ 

In (3.1) an asymptotic expansion is given that is uniform with respect to arg(z) 
for \z\ — > oo in C s . However, we observe that the proof just completed shows more: 
Allowing the reference point x to vary, the asymptotic expansion given in (3.1) is 
also uniform in x for x in a compact subset of R. 

Theorem 3.9. Assume Hypothesis 2.1, let z e C +; x e M., and denote by C £ C C + 
the open sector with vertex at zero, symmetry axis along the positive imaginary 
axis, and opening angle e, with < e < ir/2. Let M + (z,x), x > xq, be as in 
Definition 2.7. Then 

M+(z,x) = iz 1 l 2 I m + o(\z\ 1 ' 2 ) (3.58) 

uniformly with respect to arg(z) for \z\ — » oo in C s and uniformly in x as long as 
x varies in compact subsets of [xq, oo). 

Proof. Though stated for elements of the Weyl disk D + (z,xq), the proof of Theo- 
rem 3.1 shows that the asymptotic expansion given in (3.1) holds, uniformly with 
respect to arg (z) for \z\ — > oo in C e , for all elements of the Riccati disk D n {z, c, xo) 
as noted in (3.57). 

Note that the system (3.37) is independent of the reference point x . Recall 
that (5, defined in (3.34), is determined by our apriori choice of the sequence z n 
subject only to z n being in C e (c.f. (3.33)). Next note that fo(S), which is defined 
as a limit in (3.35), which is described explicitly in Corollary 3.8, and which gives 
solutions of (3.37) which satisfy (3.38) for < t < oo, is also independent of the 
reference point xq. Thus, had we chosen x' ^ xq as our reference point at the 
start, the asymptotic analysis begun in Theorem 3.7 and continued in (3.54)-(3.57) 
would remain the same after the variable change in (3.28) except for the integral 
expression present in (3.43) in which xq would be replaced by x' . However, given 
the local intcgrability assumption on Q{x) present in Hypothesis 2.1, we see that 
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the integral expression in (3.43) is uniformly continuous for x in a compact subset 
of R. Thus (3.44) and consequently (3.48) are uniform for t, and for xq, in compact 
subsets of R. 

As a consequence, we see that (3.57) holds for elements of the Riccati disk 
D^(z,c,xo), that this asymptotic expansion is uniform with respect to arg (z) for 
\z\ — » oo in C £ and that it is uniform in xq as long as xq varies in compact subsets of 
R. The asymptotic expansion described in (3.58) then follows by Theorem 3.6 □ 

We will see in the next section that the remainder term odzl 1 / 2 ) in (3.58) can 
be improved to o(l). 

4. Higher-Order Asymptotic Expansions 

In this section we shall prove our principal result, the asymptotic high-energy 
expansion of M + (z,x) to arbitrarily high orders in sectors of the type C e C C + as 
defined in Theorem 3.1. 

Throughout this section we choose Irn^ 1 / 2 ) > for z G C+. We also recall 
the following notion: x G [a,b) (resp., x G (a, b}) is called a right (resp., left) 
Lebesgue point of an element q G i 1 ((a, b)), a < 6, if dx' \q{x + x') — q{x)\ = o(e) 
(resp., /J dx' \q(x — x') — q{x)\ = o(e)) as e j 0. Similarly, x G (a, b) is called a 
Lebesgue point of q G L 1 ((a,6)) if Jl e dx' \q(x + x') — q{x)\ = o{e) as e [ 0. The 
set of all such points is then denoted by the right (resp., left) Lebesgue set of q 
on [a, b] in the former case and simply the Lebesgue set of q on [a, b] in the latter 
case. The analogous notions are applied to to x m matrices Q <E L 1 ((a,b)) mxm by 
simultaneously considering all m 2 entries of Q. The right (resp., left) Lebesgue set 
of Q on [a, b] is then simply the intersection of the right (resp., left) Lebesgue sets 
of Qj,k for all 1 < j, k < to, and similarly for the Lebesgue set of Q, etc. 

We start with the simple case where Q has compact support and hence the 
underlying half-line matrix-valued Schrodinger operator is in the limit point case 
at +oo. 

Lemma 4.1. Fix x 6 R and let x > x . In addition to Hypothesis 2.1 suppose 
that Q has compact support in [xo,oo) 7 that Q^ N ~^ £ L 1 ([x a ,oo)) mxm for some 
N G N, and that x is a right Lebesgue point of Q^ N ~ X \ Denote by M+(z,x) the 
Weyl-Titchmarsh matrix associated with Q. Then, as \z\ — > oo in C £ , M + (z,x) has 
an asymptotic expansion of the form (Im^ 1 / 2 ) > 0, z G C+) 

N 

M+(z,x) = iI m z 1 / 2 + y / m + . k (x)z- k / 2 +o(\z\- N / 2 ) 7 NeN. (4.1) 

|z|-»oo f— • 

zee. k=1 

The expansion (4.1) is uniform with respect to arg (z) for \z\ — > oo in C £ and 
uniform in x as long as x varies in compact subintervals of [x ,oo) intersected 
with the right Lebesgue set of Q^^ 1 ^. The expansion coefficients fh +t k(x) can be 
recursively computed from 

rh+,\{x) = 7^Q( X )> ™+A x ) = \Q'( x )i 
m+.k+i(x) = -(m' +k (x) + ^m +A {x)m + ^t(x)\ 1 k > 2. (4.2) 
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If one merely assumes Q G L 1 ([xo, oo)) mxm with compact support in [xo, oo), then 
(4.1) Zio^ds with N = (interpreting X\=i ' = 0); uniformly with respect to arg(z) 
/or |z| — > oo m C e and uniformly in x as long as x varies in compact subsets of 
[x ,oo). 

Proof. In the following let z e C+, Im(^ 1 / 2 ) > 0, and x > xq. The existence of 
an expansion of the type (4.1) is shown as follows. First one considers a matrix 
Volterra integral equation of the type (cf. [1, Ch. I], [24], [74], [84], [88]) 

u+(z, x) = exp(^ 1 /2 a . )7m _ J™ dx > Sm ( zl/2 J* 2 - X 'V Q( x ') u+ (z, x'), (4.3) 

and observes that the solution u + (z,x) of (4.3) satisfies u+(z, •) e L 2 ([xq, oo)) mxm 
in accordance with (3.3) and (3.4). Next, introducing 

v+(z,x) = cxp(— iz 1 ^ 2 x)u + (z, x), (4.4) 

one rewrites (4.3) in the form 

v + (z,x)=I m -^- TJ 2 J dx'[l-CK V (-2iz 1/2 {x-x'))]Q{x')v: + {z,x'), (4.5) 

and thus infers, 

M + (z, x) — u' + (z, x)u + (z, x)^ 1 = iz x l 2 l m + v' + (z, x)v + (z, x)^ 1 . (4.6) 
Iterating (4.5) then yields 

\\i + ( Zl x)\\ <C, ze C+, Imiz 1 / 2 ) > 0, x > x a (4.7) 

for some C > depending on Q. Finally, we need one more ingredient, proven in [81, 
Lemma 3] using appropriate maximal functions. Let q e £ 1 ([x ,oo)), supp(g) C 
[xo,R], for some R > 0, and suppose x £ [xq,R] is a right Lebesgue point of q. 
Then 

dx' exp(2»z 1 /2 (x ' _ x))q{xl) _JW + (\z\-^ 2 ). (4.8) 

z|^oc LIZ 1 
z€C c 

An alternative proof of (4.8) follows from [87, Theorem 1.13], which implies 

f R 

lim z- 1/2 / dx' exp(2iz 1/2 (a;' - x))\q{x') - g(a;)| = (4.9) 

|z| — >-oo J x 

zec E 

for any right Lebesgue point x of q. Given these facts, one iterates (4.5) and its x- 
derivative, integrates by parts, applies (4.8) to q = ||Qj,fe|| for all 1 < j, k < to, 
and estimates ||t5+(^, a;')|| by (4.7). Inserting the expansions for v' + (z,x) and 
v+(z, x) _1 into (4.6) (using a geometric series expansion for v+(z, x)^ 1 ) then yields 
the existence of an expansion of the type (4.1). The actual expansion coefficients 
and the associated recursion relation (4.2) then follow upon inserting expansion 
(4.1) into the Riccati-type equation (3.5). The assertion following (4.2) is an im- 
mediate consequence of (4.5) and its derivative with respect to x, (4.6), and the 
Ricmann-Lebesgue lemma. □ 

The corresponding scalar case to = 1 (for a sufficiently smooth coefficient q) 
was treated in this manner in [36]. The original version of our preprint assumed 
Q( N ) £ L 1 ([xo, oo)) mxm in Lemma 4.1. Prompted by a recent preprint by A. Ry- 
bkin [81], who used an m-function-type approach to trace formulas for scalar 
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Schrodinger operators and succeeded in removing any smoothness hypotheses on 
q, we reconsidered our original approach and extended the asymptotic expansion 
(4.1) to the case Q^ 1 ) G L^R) and x a right Lebesgue point of Q^ 1 ). 

Next we recall an elementary result on finite-dimensional evolution equations 
essentially taken from [73]. 

Lemma 4.2. ([73].) Let A G i 1 1 oc (R) mxm . Then any mxm matrix-valued solution 
X(x) of 

X'{x) = A(x)X(x) + X{x)A{x) for a.e. x G R, (4.10) 

is of the type 

X{x) = Y{x)CZ{x), (4.11) 

where C is a constant mxm matrix and Y(x) is a fundamental system of solutions 
of 

V'(x)=A(x)*(x) (4.12) 

and Z(x) is a fundamental system of solutions of 

&(x) = &(x)A(x). (4.13) 

Proof. Clearly (4.11) satisfies (4.10) since 

X' = Y'CZ + YCZ' = AYCZ + YCZA = AX + XA. (4.14) 

Conversely, let X be a solution of (4.10) and Y a fundamental matrix of solutions 
of (4.12). Define 

K(x) = Y{x)- 1 X{x), that is, X(x) = Y(x)K(x). (4.15) 

Then 

X' = Y'K + YK' = AYK + YK' = AX + YK' = AX + XA (4.16) 

implies 

YK' = XA, K' = Y- X XA = KA. (4.17) 
Thus, there exists a constant mxm matrix C (possibly singular), such that 

K(x) = CZ(x), (4.18) 
with Z a fundamental matrix of solutions of (4.13). Hence, X — YCZ. □ 

The next result provides the proper extension of Proposition 2.1 in [43] to the 
matrix- valued case. 

Lemma 4.3. Suppose Qi,Q 2 G J L 1 1 oc (K) mxm with Qx{x) = Q 2 (x) for a.e. x G 
[x ,a;i], xo < X\. Denote by Mj i+ (z,x) the Weyl-Titchmarsh matrices correspond- 
ing to Qj, j = 1, 2. Then 

(Ml + (z,x)-M^ + (z,x)) 

= -{l/2){M li+ {z,x) + M 2i+ (z,x)){M li+ (z,x) - M 2i+ {z,x)) (4.19) 
-(l/2)(M lt+ (z,x)-M 2t+ (z,x))(M lt+ (z,x) + M 2i+ (z,x)) 
for a.e. x G [x ,a;i]. 

Proof. This is obvious from (3.5). □ 
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Lemma 4.4. Suppose Qi,Q 2 S L 1 1 oc (IR) mxm and denote by Mj i+ (z,x) the corre- 
sponding Weyl-Titchmarsh matrices associated with Qj, j = 1,2. Define 

A + {z,x) = -(l/2)(M h+ (z,x) + M 2t+ {z,x)) (4.20) 

for x e [x ,xi] ; xo < Xi and assume Y + (z,x) and Z + (z,x) to be fundamental 
matrix solutions of 

$>'(z,x) = A + (z,x)ty(z,x) and $'(z,x) = $(z,x)A+(z,x) (4.21) 

respectively, with 

Y + (z,x 1 ) = Z + (z,x 1 )=I m . (4.22) 

Then, as \z\ — > oo, z G C £ , 

\\Y + (z,x )\\,\\Z + (z,x )\\ <cxp(-( 2 ; 1 -x )Ini(z 1 / 2 )(l + O (l))). (4.23) 
Proof. Define B + (z,x) by 

B+(z, x) = A+(z, x) + iz 1/2 I m , (4.24) 

then 

f Xl dx\\B + ( Z ,x)\\ = o(\z\V 2 ) (4.25) 

due to the uniform nature of the asymptotic expansion (3.58) for x varying in 
compact intervals. Next, introduce 

E+(z, x, x\) — cxp(i(xi — x)z 1 / 2 )I m , x < xi, (4.26) 

then 

Y+(z,x) = E + (z,x,xi) - j dx' E + (z,x,x')B + (z,x')Y + (z,x'), (4.27) 

J X 

Z+(z,x)=E+(z,x,x 1 )- dx' Z + (z,x')B + (z,x')E + (z,x,x'). (4.28) 

J x 

Using 

\\E + (z,x , Xl )\\ <exp(-(x 1 -x )Im(z 1 / 2 )), (4.29) 
a standard Volterra-type iteration argument in (4.27), (4.28) then yields 

|jr+(z,x )||, ||Z + (z,x )|| < cxp (-( Xl -x )lm(z^ 2 ) + J da;||B + (*,a;)||^ , 

(4.30) 

and hence (4.23). □ 

Theorem 4.5. Suppose Q 1 ,Q 2 e L j L oc (M) mxm with Qi(x) = Q 2 (x) for a.e. x € 
[x ,xi], x < x\ and let Mj i+ (z,x) be the corresponding Weyl-Titchmarsh matrices 
associated with Qj, j = 1,2. Then, as \z\ — > oo in C e , 

\\M lt+ (z, x ) - M 2 .+(z,x )|| < C(l + Izl 1 / 2 ) cxp(-2(x 1 - x )Im(z 1 / 2 )(l + o(l))) 

(4.31) 

for some constant C > 0. 
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Proof. Define for x E [x ,xi], z E C\K 

X+{z, x) = M h+ (z, x) - M 2 .+(z, x), (4.32) 

A+(z,x) = -(l/2)(M 1 ,+ (z,x) + M 2 ,+ (z,x)). (4.33) 

By Lemma 4.3, 

X' + = A+X++X+A+ (4.34) 

and hence by Lemma 4.2, 

X+(z,x) = Y + (z,x)X + (z,x 1 )Z + (z,x), (4.35) 
where Y + {z,x) and Z + (z,x) are fundamental matrices of 

x) = A + (z, x)V(z, x) and $'(2, x) = <J>(z, x)A + (z, x), (4.36) 

respectively, with 

Y+(z,x 1 )=I m , Z+{z,x 1 )=I m . (4.37) 

By Lemma 4.4, 

\\Y + (z,xo)\\,\\Z + (z,x )\\ <exp(-( a;i -a ; o)Im(^ 1 /2)( 1 + ( 1 ))) (4 . 38 ) 
as z — > exo, z e C £ . Thus, as |z| — > 00, z E C s , 

\\X + (z,x )\\ < ||X + (^a:i)|| \\Y + (z,x )\\ \\Z + (z,x )\\ 

< C(l + |z| 1/2 ) exp(-2(x 1 - x )Im(z 1 /2) (1 + o(1) )) (4.39) 

for some constant C > by (3.58), (4.35), and (4.38). □ 

The result (4.31) can be slightly improved as will be discussed in Remark 4.7. 

Given these preparations we can now drop the compact support assumption on 
Q in Lemma 4.1 and hence arrive at the principal result of this paper. 

Theorem 4.6. Fix x E K. In addition to Hypothesis 2.1 suppose that for some 
iVeN, Q^^ 1 ) E L 1 ([x ,c]) mxm for all c > x , and that x is a right Lebesgue 
point of Q( N ~^. Let M+(z,xo) be either the unique limit point or a point of the 
limit disk D + (z,xo) associated with (2.2). Then, as \z\ -tooin C £ , M + (z,xo) has 
an asymptotic expansion of the form (Irn^ 1 / 2 ) > 0, z E C+) 

N 

M+(z,x ) = iI m z 1 ^ + y / m +tk {xo)z- k / 2 +o(\z\- N / 2 ), NeN. (4.40) 

|z|-> 00 f— ' 

zec c 

The expansion (4.40) is uniform with respect to arg (z) for \z\ — > 00 in C e . The 
expansion coefficients m + ^(xo) can be recursively computed from (4.2) (replacing 

m+,k( x ) h V m +A x ))- 

If one merely assumes Hypothesis 2.1, then (4.40) holds with N = (interpreting 
Y^t=i ' — 0), uniformly with respect to arg(z) for \z\ — * 00 in C e . 

Proof. Define 




Q(x) for x E [x , xi], x < X\ 
otherwise 



and apply Theorem 4.5 with Q\ = Q, Q2 = Q- Then (in obvious notation) 

\\M+(z, x ) - M+(z, x )\\ < C(l + Izl 1 / 2 ) exp(-2(a;i - a;o)Im(z 1 / 2 )(l + o(l))) 

(4.41) 
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as \z\ — > oo, z G C e , and hence the asymptotic expansion (4.1) for M+(z, x ) in 
Lemma 4.1 coincides with that of M + (z,x ) (this also applies to the case N — 
0). □ 

Remark 4.7. A comparison of (3.58) and (4.40) with N = shows that one can 
improve the remainder term od^l 1 / 2 ) in (3.58) to o(l). In particular, this also yields 
an improvement of our rough estimate (4.39) to the effect that the term (1 + l^l 1 / 2 ) 
can be replaced by 1 in (4.39) and hence in (4.31). 

In analogy to Theorem 3.9, the asymptotic expansion (4.40) extends to one 
for M + (z,x) valid uniformly with respect to x as long as x varies in compact 
subintervals of [x ,oo) intersected with the right Lebesgue set of Q^ N ~ X \ Since 
this result represents a fundamental input for deriving trace formulas in [22], we 
provide the facts in the next theorem. 

Theorem 4.8. Fix x € M and let x > x . In addition to Hypothesis 2.1 suppose 
that for some JVeN, Q^- 1 ) G L\[x ,c)) mxm for all c > x , and that x is a right 
Lebesgue point of Q^ N ~ X \ Let M+(z,x), x > xo, be either the unique limit point 
or a point of the limit disk D+(z,x) associated with (2.2). Then, as \z\ — > oo in 
C e , M + (z,x) has an asymptotic expansion of the form (Im(z 1//2 ) > 0, z G C+) 

N 

M+(z,x) = ilrnz 1 ' 2 + y^m + . k (x)z- k / 2 +o(\z\- N ' 2 ) 1 JVeN. (4.42) 

zee. k=1 

The expansion (4.42) is uniform with respect to arg(z) for \z\ — > oo in C e and 
uniform in x as long as x varies in compact subsets of [xo, oo) intersected with the 
right Lebesgue set of Q^ -1 ^. The expansion coefficients m + ,k{x) can be recursively 
computed from (4.2) (replacing rh + .k{x) by m + .k(x)). 

If one merely assumes Hypothesis 2.1, then (4.42) holds with N = {interpreting 
J2°k=i ' = 0)> uniformly with respect to arg(z) for \z\ — > oo in C £ and uniformly in 
x as long as x varies in compact subsets of [x ,oo). 

Remark 4.9. Our asymptotic results in Theorem 4.8 are not necessarily confined to 
M + (z, a; )-matrices associated with half-line Schrodinger operators H + on [xq, oo). 
In fact, introducing in addition the analogous Weyl-Titchmarsh matrix M^(z,x ) 
associated with a half- line Schrodinger operator H_ on (— oo,xo], and noticing 
that the diagonal Green's matrix G(z,x,x) = (H — z)^ 1 (x,x) of a matrix-valued 
Schrodinger operator H = —I m -^ + Q in L 2 (R) m is given by 

G{z,x,x) = (M_(z,x) - M+(z,x))-\ (4.43) 

Theorem 4.6 then yields an analogous asymptotic expansion for G(z,x,x) in C £ of 
the form 

G(zxx) im)I m z-^ + o(\zn forJV = 0, 

[ ' ' \z\~~\(i/2)ELoGk(*)z- k - 1/2 + o(\z\- N - 1 ' 2 ) forJVGN, ^ 

where 

Go(x) = I m , Gi(x) = ig(x), etc. (4.45) 

The expansion (4.44) is uniform with respect to arg (z) for \z\ — > oo in C £ and 
uniform in x G M as long as x varies in compact intervals interesected with the 
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Lebesgue set of Q^ N ~^ (if N e N). We refer to [22] for further details and ap- 
plications of this fact in connection with trace formulas and Borg-type uniqueness 
results for Q(x), x el. 
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